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Abstract 


The following parallel algorithms are investigated in this thesis: 

• A parallel algorithm for finding k-connectcd components of interval graphs in 

O(Hoglogn) time with i og ^ sn processors on a COMMON CRCW PRAM is given, 
assuming that the intervals are given in sorted order. Further, if the end points are 
integers then it can be implemented in 0(Hog*n) time with processors on a 

Priority write PRAM. Previous algorithm for finding biconnected components took 
O(logn) time with processors on EREW model. 

• Prefix computation problem for lists is generalised to trees; the generalised problem 
is solved in O(log n) time using processors on EREW PRAM. 

• With the help of the prefix computation on trees, the time processor product of the 
algorithm for computing longest common sub-sequence [37] can be improved from 
(n log 3 n) to (n log 2 5 n). 

• Optimal O(logn) time parallel algorithms are presented for permutation representa- 
tion and coloring of cographs, given its cotree representation. The previous algorithm 
for this problem requires O(logn) time with n processors on EREW PRAM model. 

• The problem of finding OR of n bits is reduced to the problem of finding the cardinality 
of maximum independent set of a permutation graph., of n + 1 nodes. This proves a 
lower bound of fl(log n) on time for finding the cardinality of a Maximum Independent 

eyduMv^ 

Set of Permutation Graph on a Concurrent Read C o ncurr e n t Write PRAM. 

The problem of computing Parity of n bits is reduced to the problem of finding the 
cardinality of maximum independent set of a permutation graphs of 2 n nodes; thus 

viii 



proving a lower bound of ft( iogiogw ) (for the latter problem) on time on a Concurrent 
Read Concurrent Write PRAM, if only n 0 ^ 1 ) processors are to be used. 
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Chapter 1 


Introduction 


Graph theory was introduced by Euler, one of the greatest mathematicians of all time, to 
solve Konigsberg bridge problem. Since then graph theory has been used in innumerable 
number of fields — its rediscovery, by KirchhofF for analysing electrical networks and by 
Cayley for enumeration of organic chemical isomers, are some examples. A wide variety of 
applications of graphs are discussed in [20]. 

1.1 Graphs 

There is no standard terminology in graph theory. For example in [27] the author calls 
points for vertices of the graph and lines for the edges of the graph. Hence the graph 
theoretical terminology used in this thesis is first explained. 

• Graph: Graph 1 G consists of a finite set of vertices and a finite set of edges connecting 
the vertices. The graph itself is denoted by G(V, E ) where V denotes the set of vertices 
and E denotes the set of edges. If edge e,- connects node u and node v then the edge 
e,- is represented by (tt, v) or (v,u). 

• Complete graph: A graph with ^ edges, i.e. one in which every vertex is con- 
nected to every other vertex is called a complete graph. A complete graph of » vertices 
is represented by K n . 

throughout this thesis a ‘graph’ means a simple undirected graph. 
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• Vertex induced subgraph: Given a subset A C V of the vertices then the subgraph of 
G(V,E) induced by A is defined as the graph II (A, E') where E' C E and contains 
only edges from Ay. A. Tt is represented by (7([/t]). 

• Edge induced subgraph: Given a subset F C E of the edges then the subgraph of 
G(V,E) induced by F is defined as the graph 'll (A', F) where A' is the set of end 
vertices of edges in F ; i.e. v 6 A' iff there exists a node u such that («, v) € F. 

• Clique: A node induced subgraph, G([A]) of G(V,E) is a clique if and only if G([A]) 
is isomorphic 2 to the complete graph with | A | vertices. 

• Maximum clique: A clique of maximum cardinality is called maximum clique. The 
cardinality of a maximum clique is denoted by u’(G); it is called the clique number of 
G. 

• Clique cover of size k: A clique cover of size k is a partition 3 V = A x U A% U . . . U A* 
such that each A,- is a clique. The size of a smallest clique cover is denoted by k(G); 
it is called the clique cover number of G. 

• Independent set: A subset A C V of the vertices is called an independent set if and 
only if the subgraph induced by vertices of A is a null 4 graph. The cardinality of a 
maximum independent set is denoted by <*((?); it is also called the stability number of 
G. 

• Coloring: Coloring is a process of giving colors (integers) to vertices of the graph. 

• Proper coloring: A graph is said to be properly colored if no two adjacent vertices 
have the same color. Minimum number of colors needed to properly color the graph 
is denoted by x(G)] it is called the chromatic number of the graph. 

2 Graphs Gi and Gi are isomorphic to each other iff there is a one-one correspondence between their 
rertices with adjacency property. 

3 Ai nAj = </> for « jt j and Ai U . . . U Ak = V. 

4 A graph with no edges is called null graph. 
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1.2 Perfect Graphs 

The notion of a PERFECT GRAPH was introduced by Claude Berge in the early part of 
the 60’s. We know that for any graph G the minimum number of colors needed to properly 
color the vertices of G is atleast the size of the largest complete subgraph of G\ the smallest 
possible clique cover of G is atleast the size of the maximum independent set of G. More 
formally, 


<4G) < x(G) (1.1) 

a(G) < k(G ) (1.2) 

A graph G is called perfect if it satisfies the following properties: 

VACV: U (G([A])) = x(G(M)) (1.3) 

MCV:.(SP])) = *(G([A])) (1.4) 


(Note: As a(G) = iu(C?) and k (G) = x(G), if G is a perfect graph then so is G. ). 
Claude Berge conjectured that Equation 1.3 and Equation 1.4 are equivalent. Interest- 
ingly, Lovasz [36] proved the conjecture and gave the following equivalent characterization 
[35]: 


fora// A C V : u>(G([A]))a(G([A])) >| A | (1.5) 

It was customary to call a graph a-perfect if it satisfies Equation 1.4; x-P er f ec t if it 
satisfies Equation 1.3. There is no simple test for recognizing perfect graphs(see Appendix A, 
but many different classes of perfect graphs are known. Some of them are described here 
with their applications. 

1.2.1 Interval Graphs 

DEFINITION: We associate a graph G=(V,E) with the family of intervals on a real line 
!={/,■= [a;, 6,]|a,’ < 6,, 1 < i < n}, such that there is a vertex vj for each interval Ij. Two 
vertices v, and Vj are adjacent in G iff the corresponding intervals J; and Ij overlap. 
Figure 1.1(b) is an interval graph for the intervals in Figure 1.1(a). 



(a) Intervals on a real line 



(b) Interval graph corresponding to the intervals in (a) 


Figure 1.1: Interval graph 
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Figure 1.2: Channel routing 

We next look at some applications of interval graphs. 

Application 1 (Scheduling): Consider a set of courses being offered by a university. Let the 
courses be ci,C 2 , .. .,c*. Let t ,• be the time interval for the course c,-. University wants to 
assign courses to classrooms so that no two courses are conducted in the same room at the 
same time. Construct a graph G(V,E) where V is the set of courses and an edge (c,-, Cj) 
is in E if course c;’s time interval t, and course Cj’s time interval tj overlap. The graph 
constructed so is an interval graph as the time intervals are on real line. Obviously the 
problem can be solved by coloring the vertices of the graph. 

Application 2 ( Channel routing in VLSI): Components of an electrical network are laid out 
in a straight line. Certain pairs of components are to be connected using only two vertical 
runs and one horizontal run of wires see Figure 1.2. The horizontal and vertical runs are 
physically located in different layers. Each horizontal wire lies in a channel; No channel 
can carry more .than one wire at a time. The problem here is to use minimum number of 
channels. Consider an interval graph where each interval corresponds to a horizontal run 
of wire that connects a pair of electrical components. Obviously the optimal coloring of the 
interval graph gives the channel assignment using minimum number of channels. 

Interval graphs have also found applications in archaeology, biology, psychology, man- 
agement, etc. (see [25] for more details). 

1.2.2 Permutation Graphs 

DEFINITION: Let it = [iri,7r2, . . . ,x n ] be a permutation of length n then 7rf l is the 
position in the permutation where the number i can be found. Permutation graph of it, 
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Figure 1.3: Permutation graph of G([4,3,l,6,5,2]) 

G[ir] = G(V,E ) has V = {l,2,*",n} as the set of vertices and E = {(*, j)|(* — j’Xt," 1 — 
7TJ* 1 ) < 0} as the set of edges. Informally, each 7 r; is connected with 7r l+ fc if 7r,- > TTj+jt, 
for 1 < k < n — i. An undirected graph G is called Permutation Graph if there exists a 
permutation tt such that G is isomorphic to <j[ 7 t]. 

Example For the permutation t = [4,3, 1,6, 5,2] 7^ = 4, 7r2 = 3, ...7r6 = 2 and tt^" 1 = 
3, nj 1 = 6, ...TTg 1 = 4. The permutation graph G([4,3, 1,6, 5, 2]) will have vertex set 
V = {1,2, 3, 4, 5, 6} and edge set E = {(1 , 3), (1 , 4), (2, 3), (2, 4), (2, 5)(2, 6), (3, 4), (5, 6)> 
(see Figure 1.3). 

Some applications of permutation graphs are next considered. 

Application 1 (Airline): Let A and B be two sets of cities. Suppose there are air line routes 
connecting various cities in A with various cities in B, all scheduled are to be utilized at 
the same time. The problem is to assign altitudes to each flight path so that intersecting 
routes will be at different altitudes, thereby ensuring no mid air collision. Each route will 
be represented by a node in the graph G(V, E). Two nodes in the graph are connected if 
and only if the corresponding routes intersect. Coloring the nodes of the graph will solve 
our problem. It is clear (see Figure 1.4) that the graph so constructed will be a permutation 
graph. 

Application 2 (Multiprogramming): Let I={/»=[a;,6,-]|a,’ < 1 < i < A:}, be a collection 

of intervals on a line. Let | Jj | be the length of the interval /,-, i.e. 6; — a,-. Assume that 




(b) 7r-diagram 



(c) Permutation Graph 


Figure 1.4: Constructing permutation graph from airline routes 
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the intervals (they may overlap) have been sorted such that a\ < < • • • < a*. Let c,- 

be the cost associated with shifting the interval (Note: The cost is not related to the 
distance it is shifted but may depend on the length of the interval). The problem asks for 
the minimum cost shifting of the intervals such that the following conditions are preserved: 

• The order of the intervals is preserved. 

• No overlap remains, i.e. a,+i > f>,-, for 1 < i < k — 1. 

The intervals correspond to the memory requirements of k programs at a certain time in a 
multiprogramming environment; length of an interval corresponds to its cost. 

The above problem can be modeled as a graph. Each interval represents a vertex of 
the graph. There is a directed edge from interval /, to interval Ij if and only if following 
conditions are satisfied: 

• i < j 

In other words two intervals are related in the graph iff intervals between them can be 
shifted in such a way that none of these intersect. Find the maximum weighted directed 
path (fix the intervals on this path and shift the others). It can be shown that the graph 
is transitive [22]. In other words we have to find the maximum weighted clique of the 
undirected version of the above graph. Moreover the undirected version is a permutation 
graph. 

1.2.3 Cographs 

DEFINITION: A cograph is defined as follows: 

a) Single vertex graph is a cograph 

b) If Gi, G 2 , . . . , G n are cographs then their union G\ U Gi U . . . U G n is a cograph 

c) If G is a cograph then its complement G is also a cograph 

APPLICATIONS: Many graph problems are NP-Complete on general graphs. Heuris- 
tics can be used to solve these problems on general graphs. If the given input is nearly a 
special graph, but not exactly a special graph then the problem is solved with the help of 
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the algorithm for special graphs if the algorithm on special graphs is robust enough to allow 
adaptation to these near misses. 

Examination Scheduling : In this problem each course is represented by a node and node 
i, node j are connected by an edge iff some students take courses corresponding to the nodes 
i and j. Examination for the courses in the same color class can be conducted concurrently. 
In practice the graph given to this problem is close to a cograph [18]. 

1.3 Overview of the Thesis 

Chapter 2 describes the parallel computational model used in the thesis — the well-known 
Parallel Random Access Machine (PRAM) model. All algorithms in this thesis assume this 
model. Problem of finding fc-connected components of an interval graph is described in 
Chapter 3. In Chapter 4 an optimal parallel algorithm for prefix computation on trees is 
described. The algorithm takes O(logn) time with processors on an Exclusive Read 
and Exclusive Write (EREW) PRAM model. In Chapter 5 the algorithm given in [37], 
for computing a longest common sub-sequence is improved with the help of the prefix 
computation on trees 4. In the same chapter lower bounds for finding the size of the 
maximum independent set of a permutation graph are considered. As permutation graph 
is a perfect graph, the same lower bounds apply to chromatic number as well. Optimal 
parallel algorithms for permutation representation and minimum coloring of cographs are 
discussed in Chapter 6. Finally, in chapter 7 some conclusions are offered and scope for 
future work is discussed. 



Chapter 2 


Parallel Computational Model 


Parallel Random Access Machine (PRAM) model is perhaps the most popular model of 
parallel computers for algorithmic design. It neglects any hardware constraints and is 
ideal for studying inherent parallelism present in a problem. It gives absolute freedom 
to algorithm designer in presentation of parallel algorithms. In any realization of PRAM 
there will be a link between each processor and each memory location. This is however not 
realizable with the present day architectures. There are methods of simulating (like sorting 
on address and processor index) such an idealized computer on more reasonable parallel 
computers (like fixed networks of processors with number of linkages from any processor 
being bounded). This simulation costs only polylogarithmic time on most models [5]. 

PRAM is a shared memory model. Many processors work synchronously and communi- 
cate through the common random access memory. Each processor is a uniform cost Random 
Access Machine (RAM), with usual operations and instructions (see e.g., [3] for details of 
RAM model). 

The processors are indexed by consecutive integers usually from the number 0 (i.e. 0- 
(n — 1) if there are n processors). See the Figure 2.1 All processors execute the same 
instruction but (possibly) on different data. Hence it is a single instruction multiple data 
(SIMD) stream model. 

In one step each processor either reads or writes into a memory location. There are 
four models in PRAM family depending on rules for the simultaneous access of the same 
memory location. 
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Figure 2.1: PRAM model 

EREW PRAM: EREW is an acronym for Exclusive Read and Exclusive Write. In this 
model no two processors can attempt to read the same memory location at the same 
time; nor two processors can attempt to write into the same memory location at the 
same time. This is the weakest model among all PRAM models. 

CREW PRAM: This stands for Concurrent Read and Exclusive Write. In this model 
two or more processors can read the same memory location at same time. But still 
no two processors can try to write into the same memory location simultaneously. 

CRCW PRAM: In this model, simultaneous concurrent reads and concurrent writes are 
permitted. There are many variants of CRCW PRAM based on the write conflict 
resolution rule. Some of these are: 

PRIORITY: If the highest numbered (highest priority) processor succeeds in writing 
then the model is called PRIORITY concurrent write PRAM. 

ARBITRARY: If it is not known which processor succeeds, but some single pro- 
cessor is guaranteed to succeed then the model is called ARBITRARY PRAM 
model. 

TOLERANT: The model is called TOLERANT if no processor succeeds in case of 
multiple writes and contents of the memory location are guaranteed to remain 
unchanged. 
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COMMON: In this model every processor is required to write the same value in 
case of multiple writes into the same memory location. This value appears that 
location. 


See [23, 43] for simulation of one model on the other. 



Chapter 3 


Connectivity Problems for 
Interval Graphs 

We associate a graph G=(V,E) with the family of intervals J={/,=[a,-, 6,- ]|a; < 6j, 1 < i < n}, 
such that there is a vertex Vj for each interval Ij. Two vertices v,- and Vj are adjacent in G 
iff the corresponding intervals /,• and Ij overlap. 

Interval graphs have found applications in archaeology, biology, psychology, manage- 
ment, VLSI design, scheduling, etc. (see [25] for more details). 

Efficient parallel algorithms have been devised for various problems for interval graphs 
[11, 42, 47]. In this thesis an O(loglogn) time algorithm with ]fc --- processors for finding 
biconnected components is obtained on a COMMON CRCW PRAM, assuming that the in- 
tervals are given in sorted order. This algorithm can also be implemented on a TOLERANT 
CRCW PRAM with same resource bounds. We say /,• > Ij iff a,- > aj. The intervals are 
sorted in increasing order if I\ < I 2 < • • * < 7 n . Further, if the end points are integers then 
the time to find biconnected and k-connected components can be reduced to 1 O(log* n) and 
O(fclog*n), respectively with processors on a Priority write PRAM. 

Earlier Liang, Dhall and Lakshmivarahan [34] have obtained O(logn) time algorithm 
for finding biconnected components with n processors on a CREW PRAM. Sprague and 
Kulkarni [45] have given O(logn) time optimal (on sorted intervals) algorithm for finding 
biconnected components on an EREW PRAM. Das, Calvin and Chen [19] have obtained 
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O(logn) time algorithms for finding all bridges and articulation points with n processors 
on an EREW PRAM. 

We generalize the notion of “prefix- maxim a” to “k-th-prefix-maxima”. All algorithms 
are based on the observation that k-connected components can be easily identified if “k- 
th-prefix-maxima” are known. We believe that the notion of “k-th-prefix-maxima” will be 
useful in other applications also. 

3.1 K-th Prefix Maxima 

Given a collection of intervals 2= {[a,-, 6,] |1 < * < n). We say a; and 6,- are the left end point 
and the right end point of i th interval respectively. We will say two intervals are adjacent, 
if the corresponding vertices in associated graph are adjacent. 

Right Interval: I T is called the right interval of set S if h r — min{6;|i € <S}. 

Left Interval: Ii is called the left interval of set S if a/ = max{a,ji € «S}. 

We define k-th-prefix-maxima(Ii) to be the k th largest element in the set 5,- = {bj\aj < a,}. If 
k-th-prefix-maxima/( /,•) is 6 m then we say interval I m contributed to k-th-prefix-maxima(Ji). 
We say interval Ii is critical in the j th iteration iff j-th-prefix-maxima(Ii) is greater than 
j-th-prefix-maxima(/,_i). We say that a set of intervals I={/,} spans [a,b], if for every 
point x, a < x < b, there is an interval Ij, such that aj < x < bj 

The k-th-prefix-maxima (for k > 2) of each interval can be found as follows. 

Step Ij Let C[j]:=6j for 1 < j < n. 

Step 2: Find prefix maxima of the array C. Let array D be the prefix maxima 
of array C. 

Step 3: For each critical interval Ij (i.e. D[j]^D[j-l]) do C[j]:=D[j-l]. 

Step 4:Repeat step(2) followed by step(3), k - 2 more times. 

Step 5: Find prefix maxima of array C. 

The last prefix maxima will give us the k-th-prefixmaxima. 

Prefix maxima can be found optimally (with a processor-time product of 0(n) ) in 
O(logn) timeon an EREW PRAM, and in O(loglogn) time on a COMMON CRCW PRAM 
[9, 8]. If the items are integers then prefix maxima can he found in optimal O(log* n) time 
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on a Priority write PRAM [8, 9]. 

The usual “1-st-prefix-maxima” will be called prefix maxima. 

3.2 Connected Components 

We compute connected components using 2-nd-profix-maxima. Note that if for any interval 
2-nd-prefix-maxima(/;) < a,- then the l-st-prefix-maxima(/,) will be 6; itself. Hence there 
is no interval spanning the real line from 2-nd-prefix-maxima(/,') to a,-. Or, /,• starts a new 
connected component. 

All intervals which start a new component can be identified (2-nd-prefix-maxima(/i) < a,) 
and marked. Let A[i] be oo if interval R is not marked, and let A[i]=i, if interval U is marked. 
Then find the left nearest smaller of array A. [10, 51]. Thus each interval is identified by 
the smallest interval in that component. The above procedure takes O(loglogn) optimal 
time on a COMMON CRCW PRAM and O(logn) optimal time on an EREW PRAM. 

Theorem 3.1 Connected components of interval graph can be found in O(loglogn) optimal 
time on a COMMON CRCW PRAM. I 

We can do better on a PRIORITY write PRAM, if the end points are integers. We 
proceed as follows 

A flag is set to 1 if the interval is marked and it is reset to 0 otherwise. We 
“connect” all l’s by using ordered chaining algorithm of [41] in 0(a(n)) optimal 
time on a COMMON CRCW PRAM. Each 0 can find the index of the imme- 
diately previous 1 by using the algorithm of [9] for the nearest one complement 
problem. This operation also takes 0(a(n)) optimal time on a COMMON write 
PRAM. 

Theorem 3.2 Connected components of interval graph can be found in 0(log* n) time with 
~ processors on Priority write PRAM if the end points are integers. I 

Note that if all end points are given in sorted order, then we can replace each end point 
by an integer. 
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3.3 Biconnected Components 

Biconnected components of a connected graph G are found by eliminating all articulation 
points of G and finding connected components of the resulting graph G'. Articulation points 
are found by using 2-nd-prefix-maxima algorithm. 

If 2-nd-profixmaxima.(/j_i) < a,' then the interval which contributed to prefix maxima is 
an articulation point, provided 2-nd-prefix-maxima(J,_ 1 ) > 0. This can be seen as follows. 
Let us assume that If contributed to l-st-prefix-maxima(J,_i) and I a to the 2-nd-prefix- 
maxima(7,_i). Thus 2-nd-prefix-maxima(/,_i) = b 3 < aj. Moreover, as G is connected, 
bf > ai . Thus, the real line from b a to a; is spanned by only one interval, If. Hence, If is 
an articulation point. 

We identify and then remove all the articulation points of G, and find the connected 
components of the resulting graph G' . In more detail, we proceed as follows. 

If interval /,• is an articulation point for components Cj\ , Cj2, C/ 3 , • • • , Cj m then every 
component’s first interval (smallest numbered interval in that component) except the first 
component Cj\ detects the articulation point (using the prefix maxima routine). Then, all 
detected articulation points are removed, and connected components are found using the 
algorithm of Section 3.2. We again, “link” all first intervals of each connected component 
as in Section 3.2; each component corresponds to a bi-connected component. We next 
describe the method for recording each articulation point in all components (to which it 
is adjacent). Recall, that every component’s first interval, except the first component Cj\ 
knows the articulation point. We still have to inform the first component. 

Observe that each bi-connected component can be adjacent to at most two articulation 
points, one adjacent to first interval of the component, and other to the interval which ends 
last. We will use an array A of length 2n to record articulation points. If Ij is first interval of 
a bi-connected component, then the articulation points will be available in locations A[2j-1] 
and A[2j]. 

• If interval /,• detects articulation point Ij then record Ij at location 2 i in the array A. 

• Let Ik be the first interval of a component just before “component” /,• then record artic- 
ulation point Ij at location 2k -f 1 in the array A. 

Thus, each interval knows the bi-connected component in which it lies (first interval 
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of the bi-connected component) and each first interval Ij has pointers to at most both 
articulation points. This, completely characterizes the biconnected components. 

Theorem 3.3 Biconnected components of interval graph can be found in 0(loglogn) time 
optimally on a COMMON CRCW PRAM model. I 

Theorem 3.4 Biconnected components of interval graph can be found in 0(log* ») time 
optimally on Priority write PRAM if the end points are integers. 1 


3.4 k-Connected Components 

Given a positive integer k, a graph G=(V,E) with at least k + 1 vertices is called k-vertex 
connected if deletion of any k — 1 vertices leaves the graph connected. Testing k-connectivity 
of a general graph [13] with n vertices and m edges takes 0(fc 2 logn) time using 
processors on CRCW PRAM. 

For simplicity, we assume that the graph G is k — 1 connected. We say a 6et of vertices 
form a disconnecting set if removal of all vertices in the set disconnects the graph. As 
the graph is k — 1 connected, it follows that every pair of vertices in a disconnecting set 
are adjacent; and in each component adjacent to the disconnecting set there is at least a 
vertex which is adjacent to all vertices of the disconnecting set. Further, each k-connected 
component can be adjacent to at most two different disconnecting sets. We will find k- 
connected components by eliminating all maximal disconnecting sets of size k — 1. 

Observe that if k-th-prefix-maxima(Ii-i) < a; then all those intervals which contributed 
to l-st-prefix-rnaxima(/,-_i), 2-nd-prefix-maxima(/,-_i), • • •, k - l-th-prefix-maxima(/j_i) 
will form a disconnecting set of size k — 1. 

Lemma 3.1 For each disconnecting set of size k — 1 there exists an interval j such that 
k-th-prefix-maxima(Ij-i ) < Oj. 

Proof: Let V be a disconnecting set of size k - 1. Let I 7 be the set of intervals obtained 
by removing all the intervals of the disconnecting set V. Let C\ and Ci be the first two 
components of the resulting graph. Let I r be the right interval of C\ and let Ij be the first 
interval of (7 a . As C\ and (7a are in different components in I', prefix-maxima'(a r ) < aj. 
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As there k — 1 intervals in disconnecting set, it follows that k-th-prefix-maxima(a r ) < a,-, 
and for any Iy € V k-th-prefix-maxima(a J /) < a,j. Thus, the lemma follows. 1 

Thus, all k-connected components can be found by first removing all disconnecting sets 
of size k — 1 and then use connected components algorithm on the resulting graph; connected 
components of the resulting graph will correspond to k-connected components of the original 
graph. Vertices of disconnecting sets can be made known to the components adjacent to 
them by a method similar to that described in Section 3.3. However, as there are A: — 1 
vertices in disconnecting set, we will have to use k — 1 arrays A\,- • - Ah, each of length 2n. 

Theorem 3.5 k-connected components of interval graph can be found in 0(k log log n) time 
with ( i og i^ gw ) processors on a COMMON CRCW PRAM model. I 

Theorem 3.6 k-connected components of interval graph can be found in 0(fclog*n) time 
with ( t og * n ) processors on Priority write PRAM if the end points are integers. 1 



Chapter 4 


Prefix Computation on Trees 


Prefix computation on lists and arrays is a well known problem; it is a very useful routine 
in graph theoretical problems (see e.g., [33]). Generalization of this problem on trees is: 

Given a tree T and an associative binary operation 0, compute for each node v: 
vj 0 u 2 O • • • 0 vjt; here is parent(v), v* is the root of T and v, + i is parent(w,), 
for 1 < i < k — 1 . 

Prefix computation problem for linked list can be reduced to the prefix computation 
problem for arrays using list ranking. List ranking can be solved in O(logn) time with 
processors on an Exclusive Read Exclusive Write (EREW) PRAM [16, 6]. Prefix 
computation on an array can be done optimally in O(logn) time (see e.g., [33]) on an 
EREW PRAM. 

Remark : Another generalization of the prefix computation problem is the so called “General 
prefix computation” problem [46]: Let X , Y be two given arrays and an associative binary 
operation © be defined on elements of X ; elements of Y can be compared by a linear order 
<. Then for each index t we are required to compute Z{ = Xj t 0 Xj 3 © • ■ • © Xj k where 
Ji < h < < jk < i and {j\, j 2 , . . ., j m , . . .jk} is the set of indices for which Y jm < Yi 

for 1 < m < k. In other words for each X,- we “add” only those Xj for which a ) j < i b) 
Yj < Y{. Sorting can be reduced to general prefix computation problem (Y consists of items 
to be sorted and X consists of n l’s). An O(logn) time, n processor parallel algorithm for 
this problem has been obtained in [46]. 
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Remark : An O(logn) time parallel algorithm for Prefix Computation on Trees with n 
processors can be obtained by using the “General Prefix Computation” algorithm of [46]. 
Traverse the tree in pre-order and pack the given items into array X according to their 
pre-order number; array Y contains the rep-order number of nodes; “rep-order” traversal 
can be described as follows: visit a node, then visit the children in right to left order. Rep- 
order number can be found by first reversing the order of children of each node and then 
traversing the tree in pre-order. 

4.1 Algorithm 

The given tree T is first converted into a regular binary tree B (each node has two or no 
children) using the method of [32]: 

If a node v has k children, v \ , V 2 , . . . , u* then introduce k additional nodes v [ , v ' 2 , . . . , v' k such 
that 

• is a right child of v. 

• Vj is a right child of v '-_ x , for 2 < j < k 

• Vi is a left child of v. 

• Vj is a left child of v'_j , for 2 < j < k 

Note that the ancestors of a node in T are also ancestors of the corresponding node in 
B. But the set of ancestors of a node in B may include some new nodes introduced during 
the conversion. If the identity element of the operator ©, is put in the newly introduced 
nodes, the problem of prefix computation on T can be solved once the problem is solved for 
B. Further, note that the number of nodes introduced is always less than n. Thus, without 
loss of generality, in the rest of the chapter, we may assume that the given tree is a binary 
tree. 

The algorithm successively modifies the given tree T into a sequence of “simpler” trees 
Ti, T 2 , . . . , Tk such that 

a) the final tree T* is simple enough to obtain prefix computation in 0(1) time 

b) the prefix computation on T< can be done in 0(1) time using the “knowledge” of prefix 
computation on !T,-+i and 

c) k is small, i.e., k = O(logn). 
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All those conditions are satisfied for the tree contraction sequence generated by algorithm 
of [1]. The tree contraction sequence is defined as follows[l]: 

DEFINITION : Let B be a binary tree with vertex set V(B). A sequence of trees 
B\ , B%, • • • , Bk is said to be a tree contraction sequence of length k for B if, 

1. J?i = B 

2. V(Bi)CV(Bi. x ) 

3. | V(B k ) |< 3 

4. 1ft; € - V(Bi) then either 

i) w is leaf of fi;_i or 

ii) v has exactly one child, x, in x € V(B{), and the parent of v in 

Bi - 1 is the parent of x in B{. 

We use the basic algorithm of [1] for tree contraction. How-ever some extra bookkeeping 
information is maintained in routines — prunef ] and bypass[ ]. 

Algorithm: We assume that the field “Valuefv]” stores the value initially at node v, or the 
identity element of “O”. Number the leaf nodes from left to right s.t., for left-most leaf v/, 
Number[u/]=0 and for right-most node v T , N umber [u r ]=n - 1. The algorithm is as follows: 

for i=l to flogn] — 1 do serially 

{ 

for each leaf v pardo 
{ u;:=parent(t;); 

if Number(t;) is odd and w ^ root 

if v is a left child of w then j:=0; prune(v); j:=l; bypass(tt;); 
if v is a right child of w then j:=2; prune(v); j:=3; bypass(tu); 
else Number(v):=Numbcr(v)/2; 

> 

} 

JENTk'AL 1:17 A RV 
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prune(v) 

{ 

Index[v]:= 4 i + j; 

Last_Parent[u]:= Parent[i>] 

Remove the node v as in the usual “Prune” routine 

} 

bypass(u) 

{ 

Index[u]:= 4 i + j; 

Last_Parent[u]:= Parent[v] 

If x is the only child of v then 
parent(x) := parent(u) 
value(x) := value(x)© value(u) 

Remove the node v as in the usual “Bypass” routine 

} 

Note : If Index[v] is “/” , then the node v was in trees Bj, Bi, but was not in tree B/ +1 . 

Lemma 4.1 Let f? t be the i th tree in the tree contraction sequence, w be the parent of node 
v in Bi and vi,v 2 , .. .,Vj be the proper ancestors of v and proper descendants of w (i.e. 
v ii w 2 i • • * > Vj are the nodes on the path from v to w) in B. Before computation of B,+ j, 
the value (v) will be Init-value(v) © Init.value(vi) ©•••© Init.value(vj) where Init.value 
represents the initial value given to the nodes. 

Proof: : The lemma is proved by induction. The lemma is trivially true for B\ . Assume 
that it is true for f? t . 

Assume that node u is bypassed in B, to get B, +1 . Let v be the only child of u and w 
be parent of u in B,-. By induction hypothesis 
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prune(u) 

{ 

Index[v]:= 4 i + j\ 

Last -Parent [v]:= Parent [v] 

Remove the node v as in the usual “Prune” routine 

} 

bypass(v) 

{ 

Index[u]:= 4i 4- j; 

LastJParent[u]:= Parent [u] 

If x is the only child of v then 
parent(x) := parent(t;) 
value(x) := value(x)© value(v) 

Remove the node v as in the usual “Bypass” routine 

} 


Note : If Index[u] is then the node v was in trees J?i, . . . , J3/, but was not in tree -Bz+i- 


Lemma 4.1 Let Bi be the i th tree in the tree contraction sequence, w be the parent of node 
v in Bi and vi , v?, ■ ■ - ,vj be the proper ancestors of v and proper descendants of w (i.e. 
Vi,V 2 , . . . , vj are the nodes on the path from v to w) in B. Before computation of J9;+ 1 , 
the value(v) will be Init.value(v) © Init. value (vi) ©■•■© IniLvalue(vj) where IniLvalue 
represents the initial value given to the nodes. 

Proof: : The lemma is proved by induction. The lemma is trivially true for B\. Assume 
that it is true for i?;. 

Assume that node u is bypassed in Bi to get . Let v be the only child of u and w 
be parent of u in J?,\ By induction hypothesis 
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value(t>)= Init_value(w) © Init-value(vi) ©•••© Init_value(u m ) and 
value(u)= Init_value(u) © Init_value(ui) ©•••© Init_value(u n ) 

where nodes vi , w 2 , . . . , v m are on the path from v to u and , i> 2 , • • • , u n are on the path 

from u to w. After bypassing node u (see bypass routine), value(v) will be 
value(u)= Init_value(u) © Init_value(ui) ©•••© Init-valuefum) © Init_value(u) 0 
Init-value(ui) ©• ••© Init_value(u n ). 

All these nodes are on the path from v to w in B ( w is parent of v in I?,+ 1 ). I 

Corollary 4.1 Let Index[v] be i and w=LasLparent[v], be the parent(v) in Bi. When the 
algorithm terminates the value(v) will be 
Init-value(v) 0 Init.valuef v\) ©•••© IniLvalue(vj) 

where vi,v %, . . . , Vj are nodes on the path from v to w in B. 1 

The last tree in the sequence, Bk has three nodes, r, the root of tree, ci, left child of r and 
C 2 , right child of r. By Corollary 4.1, final value of ci and C 2 can be calculated as follows: 
value(cj) := value(cj)© value(r) 
valuc(c 2 ) := value^)© value(r) 

Assume that the final value of each node in Bi+i has been calculated. The value of 
nodes in Bi but not in Bi+i can be calculated as follows: Let v be a node in 2?,- but not in 
Bi + i \ i.e., Index[v]=i. Let w=Last_Parent[v], be the parent of v in Bi then value(tw) will be 
Init_value(tu) © Init_value(i«i) ©•••© Init_value(u>j) where w\, W 2 , • • • > w T are proper an- 
cestors of w in B and value(v) will be 

Init_value(v) © Init_value(tJi) ©•••© Init_value(vj) where V\,V 2 3X6 nodes on the 
path from v to w in B. Final value(u) will be: value(u):=value(v) © value(tn). 

As the maximum value of “Index[ ]” is 41ogn, the algorithm can easily be implemented 
in O(logn) time with n processors. Using sorting algorithm of Cole and Vishkin [15], 
processor-allocation problem can be solved and the number of processors reduced to 
without any increase in time. In more details, we proceed as follows: As “Index[ ]” is in the 
range 1, . . . , 4 logn, the nodes can be sorted on the key Index[v], in O(logn) time with 
processors [15]. Calculate each node's RANK as, RANK[v]:= POSITION(r)+mj^j- where 
m=Index[v] and POSITION(v) is the position of node v in the sorted array. We process 
nodes in 2 logn rounds — In j tfi round all vertices whose RANK lies between (2 logn — 
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j — l)(p^) and (21ogn — are processed. This ensures that no two vertices with 

different Index[ ] (or in different trees) are processed in the same round. 

Let Ci be the number of nodes in B, but not in B,+i. Time taken in i th round is 
majc(^*^ si , 1); thus time taken by algorithm is atmost 

E-I°i gn maz(£^, 1) < + 1) = '-T ES sn Ci + 4 log n = log n + 4 log n = 

51ogn = O(logn) 

Thus, the following theorem follows: 


Theorem 4.1 Prefix computation on trees can be done in O(logn) time using pro- 
cessors on an EREW PRAM. 1 
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4.2 Applications 

Applications of prefix computation on trees are discussed in Chapter 5 and Chapter 6. 
This algorithm has applications for computing a largest increasing sub-sequence which in 
turn can be used to compute maximum independent set of a permutation graph. Previous 
algorithm for computing largest increasing sub-sequence[37] takes 0( log 3 n) time with n 
processors on a CREW PRAM. The processor-time product can be reduced from 0(n log 3 n) 
to 0(n log 2-5 n) on a CREW PRAM; see Chapter 5 for details. Prefix computation on trees 
is also used for coloring and permutation representation of cographs; see Chapter 6. 



Chapter 5 


Maximum Independent Set of a 
Permutation Graph 


5.1 Introduction 

let 7 r = [7ri,7T2,. ..,7r n ] be a permutation of length n then 7T," 1 is the position in the 
permutation where the number * can be found. Permutation graph of 7r, G[7r] = G(V, E ) 
has V = {1,2, ■ ■ -,n} as the set of vertices and E = {(*, j)|(* - - < 0} as the 

set of edges. 

Informally, each 7r; is connected with ir, + fc if 7 r; > 7r, + jt, for 1 < k < n — i. An undi- 
rected graph G is called Permutation Graph if there exists a permutation tt such that G is 
isomorphic to <3[7r]. 

Example For the permutation 7r = [4, 3, 1,6, 5, 2] t\ = 4,7r2 = 3,...7rg = 2 and 7T] -1 = 
3,7^* = 6,. . . TTg"* = 4. The permutation graph G([4,3, 1,6, 5, 2]) will have vertex set V = 
{1,2, 3, 4, 5, 6} and edge set E = {(1, 3), (1,4), (2,3), (2,4), (2, 5)(2,6), (3,4), (5,6)>. I 

Permutation graphs were introduced by Even and Pnueli in 1972 [21]; for applications of 
permutation graphs see e.g.[25, 22]. Spinrad gave 0(n 2 ) algorithm for recognizing permu- 
tation graphs [44]. Finding maximum independent set of a general graph is NP-complete 
[24], but maximum independent set of a permutation graph can be found in polynomial 
time [31, 52]. 

Next section describes the algorithm for finding maximum independent set of a permu- 
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tation graph. Lower bounds for finding maximum independent set of a permutation graph, 
given its permutation representation, are presented in Section 5.3. 

5.2 Algorithm for Maximum Independent Set 

String .S’ is said to be a sub-sequenco of siring A if either a) .S’ is equal to A or b) S is 
obtained by deleting some symbols from A. String C is a common sub-sequence of string 
A and string B if C is a sub-sequence of A as well of B. A string having maximum length 
among all common sub-sequences is called a longest common sub-sequence. The longest 
common sub-sequence problem can be used for computing maximum independent set of 
a permutation graph; for other applications of longest common sub-sequence problem see 
[38, 14]. 

5.2.1 Preliminaries 

Sub-sequence S = -Si, -S 2 , • • • , s r of a string A is called an increasing (decreasing) sub- 
sequence if s, + i > Sj (s, + i < s,) for 1 < t < r — 1. Increasing (decreasing) sub-sequence of 
string A which has maximum length among its increasing (decreasing) subsequences is called 
Longest Increasing (Decreasing) Sub-sequence. It is known that [25] increasing (decreas- 
ing) sub-sequence of a permutation corresponds to independent set ( clique) of the graph 
represented by the permutation. Hence a Longest Increasing (Decreasing) Sub-sequence 
corresponds to a Maximum Cardinality Independent set (Maximum cardinality Clique) of 
the corresponding permutation graph. 

Largest Increasing sub-sequence can be obtained by finding the longest common sub-sequence 
of the given string and the string [1,2, •••,«]. 

5.2.2 Algorithm 

The algorithm of [37] for finding the Longest Common Sub-sequence takes O(log 3 n) with 
max(r, n) processors, where r is the total number of positions at which the two input strings 
match and n is the length of the longer string. Some of the steps in the algorithm can be 
implemented at lesser cost with the help of prefix computation on trees . For completeness 
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Figure 5.1: Lu’s Algorithm 

the algorithm of [37] is breifly described below (for details refer to [37]). 

Let A = ai,a2,"-,a n and B = 61,63, •••,6 m be the two strings. Let G[i,j] be the 
length of longest common sub-sequence of strings a\ , . . . , a,- and 61 , . . . , bj. Only matched 
symbols are relevant in constructing a longest common sub-sequence. G[i,j] is computed 
only if a; = bj for 1 < * < n; 1 < j < m. The selected entries, where a,- = bj, on the grid 
(see Figure 5.1) are called points and the value G[i,j] is called class of point (i, j). Point 
q(i q ,jq) is said to be dominated by point p(i p ,j p ) iff i p < iq and j p < j q . Point f(ij,jj) 
can be the father of point s(i B ,j t ) if 

1. point / dominates point s 

2. class(/) is one less than the class(s). 

A processor is assigned to each pair ( i,j ) if a, = bj. This scheduling step can be done 
in 0(logm + logn) time using m + n processors (see [37]). 

Divide the grid into two parts by a horizontal line. Solve the problem (i.e., find the class 
of each point in that half) for each half concurrently and merge the two solutions. After 
the merge the class of each point in the upper half will remain the same, but the points in 
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the* lower half might have got their class changed. Hence, in the merge step wo have to find 
the right class for each point in the lower half. 

The merge step is carried out as follows: Project the left most point of each class in 
the upper half onto the horizontal dividing line. Denote the projection as joint. Locate 
a dummy joint at the left side of the grid on the horizontal line. Merge consists of two 
phases. In phase I each point in the lower half finds its father. A point in class 1 finds its 
father among the joints. Hence a forest is formed with joints as the roots. Roots will inform 
their class numbers to their descendants. Each point can decide its new class in phase I by 
adding its root’s class to its depth in the tree. This phase is called update with joint. 

The second phase of the merge step is update between trees i.e. merging the trees formed 
in the first phase. This phase consists of log t iterations if there are t trees. In iteration i, 
2‘ trees should be merged. Let us call the left 2 ,_1 trees as left trees and similarly the right 
2 ,_1 trees as right trees. Each point p finds a point t in the right trees such that the point 
t dominates point p and point t has maximum class number among such points. Then the 
point p changes its class number to ma x(class(p),class(t) + 1). The following is a brief 
description: 

Algorithm LCS 
Begin 

Divide the input grid into two parts, upper half and lower half 
Recursively find the class numbers of each part independently and concurrently 
Merge the two parts 
End. 


Procedure Merge 
/* Update with Joint*/ 

Begin 

Find Joints from the upper half 

Sort the points in the lower half lexicographically on the key (class number, column number) 
Each point in lower half finds its father in the class above it; 
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however a point in class 1 finds its father from the Joints. Thus forming trees. 

Perform propagation _1 (class) on each tree; i.e. each point p gets class(root(p)) + depth(p) 
/* Ujxlatc between trees*/ 

For logo- + 1 times do serially /* a — number of Joints*/ 

Each point p(t p ,j p ) finds the point q(i g ,j q ) in the right trees such that 

j p > j q and q has the maximum class number among the points that dominate p. 
diff(p)=class(q)-class(p) 

Perform propagation_2(diff); i.e. each point gets maximum diff among its ancestors 

class(p):=class(p)+difF(p)+l 

if difF(p) > diff(father(p)) then father(p):=q 

End 


First improvement comes from the fact that we can use integer sorting instead of com- 
parison based sorting. Integer sorting of n numbers in the range [1 • • • n] takes 

1. 0(logn-v/logn) with (j^) processors on an EREW PRAM [4]. 

2. 0( | 0 ~ ^n time with processors on a ARBITRARY-write CRCW PRAM 

[ 12 ]- 

Operations propagation.1 and propagation^ can be performed in O(logn) time with 
(logn) processors on an EREW PRAM as follows: 

Note that class of each point will be its level number in the tree. Class of root (Joint) 
can be propogated by assigning n + 1 to every other node and then performing prefix 
computation on the tree 4 with minima as the associative binary operation. The class of a 
node p after propagation-1 will be class(root)+class(p). Similarly, propagation.2 can again 
be done by prefix computation on tree with maxima as the associative binary operation. 

Algorithm for prefix computation on trees assumes Adjacency List Representation. We 
can sort the edges using Integer Sorting algorithm [4, 12] to convert the representation to 
adjacency list. 

Hence the following theorem, 
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Theorem 5 . 1 The algorithm for finding a longest common sub-sequence of two given strings 
takes 

1. O^og 3 \/logn) time with processors on an EREW PRAM model 

2. 0{ log 3 ) time with max( n , r ) processors on an ARBITRARY-write 
CRCW PRAM 

where r is the total number of positions at which the two given strings match. 

Corollary 5.1 The algorithm for finding a maximum cardinality independent set or a max- 
imum cardinality clique of a Permutation graph takes 

1 . 0(log 3 riy'log n) with (j^-) processors on an EREW PRAM. 

2. OClo^n) time with processors on an ARBITRARY-write CRCW PRAM. 


5.3 Lower bounds for the maximum independent set 

Given a binary string 1 (with alphabet {0, 1}) A = • •,<*„, we construct permutation 

x by looking for instances of 10’s (that is digit 1 immediately followed by digit 0) in A as 
follows: 

If a,- = 1 and a,+i = 0 then let x,- = i + 1 and x,-+i = i. In all other cases let x,- = i. 
Lemma 5.1 x is a permutation. 

Proof: We first claim that every integer is present in the sequence. For the sake of argument 
assume that integer j is not present in the sequence x. If aj = 1 and aj+x = 0 then x J+ j is 
j. If aj = 0 and aj-i = 1 then xy_i is j. In all other cases Xj is j itself. Thus every integer 
in the range [l...nj is present in the sequence. 

No integer is repeated as the length of the sequence x is only n. Hence x is a permutation. 

I 

Theorem 5.2 In 0(1), time with n processors on an EREW PRAM, the problem of finding 
OR of n bits can be reduced to the problem of “ detecting ” presence of pattern 10 in a string 
of length n + 1. 


1 In this chapter we will use the term “string” for binary string 
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Proof: Let the given input string be X = Create another string Y - 

yuV2,- ■■■.Vn+i where 



if 1 < i < n 
if * = n 4- 1 


OR of given input is TRUE iff there is atleast one 10 pattern in Y. 


I 


Theorem 5.3 In 0(1), time with 2 n processors on an EREW PRAM, problem of finding 
Parity of n bits can be reduced to the problem of finding the parity of the number of 10 
patterns in a string of length 2n. 


Proof: Let the given input string be X - x x ,xi,- Create another string Y = 

V \ , V 2 , • • • , V 2 n with p 2i _i = n and y 2 « = 0 for 1 < t < n. 

Parity of given input is EVEN if and only if the no. of 10 patterns in the string Y is EVEN 
(In fact the no. l’s in the string X is same as the number of 10 patterns in the string Y). 
Hence the result. I 


Lemma 5.2 There is a one-one relation between 10 patterns in the given string and the 
edges in the permutation graph of x constructed from the string. 

Proof: Xi can be either i - 1, t, or s + 1. Thus x,- < T,- + jt, for 3 < k < n - i. Further if 
7T; is i + 1 then x,+i = i; moreover as x is a permutation, x t+ 2 can not be i + 1 Hence, we 
can say that x,- < Xj+jt, for 2 < k < n — i. This implies that there is no path of length 2 or 
more in the permutation graph. In other words any node * is connected with atmost one 
node — node i - 1 or node i + 1 (but not both). 

Let ajOj+i be an instance of 10 (i.e. «,• = 1 and a, +1 = 0). Then x,- = * + 1 and x,+i = » 
hence there is an edge in the permutation graph between nodes i and i + 1. 

Conversely, if there is an edge between nodes j and j + 1 in the permutation graph, 
then Xj = j + 1 and Xj+i = j, or aj = 1 and aj+i = 0. Hence for every edge, there is a 10 
pattern. ® 

Theorem 5.4 In 0(1) time, problem of counting the number of 10 patterns can be reduced 
to the problem of finding the cardinality of maximum independent set of permutation graph. 
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Proof: Consider the graph of permutation x constructed from the given string. Each vertex 
in the graph is either of degree 0 or 1. Let n 0 be the number of vertices of degree 0, and 
Tii be the number of vertices of degree 1; n 0 + nj = n. Let e be the number of edges in the 
graph; e = Let I be a maximal independent set. Clearly all n 0 vertices of degree 0 will 
be in I. As no vertex is of degree 2 or more, exactly one end-point of each edge will be in 
I. If CM IS is the cardinality of maximum independent set, then CM IS = n 0 + e = n - e. 

But by Lemma 5.2 the number of 10 patterns in the graph is exactly equal to e. Thus, 
if we know the number of 10 patterns we know the cardinality of maximal independent set. 

I 

As problem of finding “OR” of n-bits requires fl(logn) time on a Concurrent Read 

Exclusive Write PRAM [17], we get the following corollary. 

■ 

Corollary 5.2 The problem of finding^ maximum independent set of a permutation graph 
of n nodes requires fl(logn) time on a Concurrent Read Exclusive Write PRAM. 

As problem of finding “PARITY” of n-bits requires fl( | o ^” n - ) on a Concurrent Read 
Concurrent Write PRAM [7], we get the following corollary. 

Corollary 5.3 The problem of finding cardinality of a maximum independent set of a per- 
mutation graph of n nodes requires time on a Concurrent Read Concurrent Write 

PRAM, if only n°^ processors are to be used. 

Remark: It is known that the complement of a permutation graph is also a permutation 
graph. In fact it can be shown that 

G[x r ] = G \ W 

where ir r is the permutation obtained by reversing the sequence x and G[x] represents 
the complement of the graph G[x]. As the cardinality of maximum independent set of G 
equals the cardinality of maximum clique of C? , the lower bounds proved above apply to the 
problem of computing the cardinality of maximum clique of permutation graph. 



Chapter 6 


Permutation Representation and 
Coloring of a Cograph 


Many graph problems like coloring, maximum independent set, minimum dominating set 
etc. are NP-Complete or intractable for general graphs [24]. Most of these problems have 
a polynomial time algorithm for special graphs like Interval graphs, Permutation graphs, 
Cographe and Chordal graphs (see [25] for definition). 

Cographs are also called Digraphs [30], P 4 restricted graphs [18] and Hereditary Dacey 
graphs [48]. Several characterizations are known for cographs, but perhaps, the most im- 
portant characterization for algorithm design is its cotree representation. 

6.1 Definitions and Preliminaries 

DEFINITION: A cograph is defined as follows: 

a) A graph consisting of a single vertex is a cograph 

b) If . . . ,(7„ are cographs then their union U G? U . . . U G n is a cograph 

c) If G is a cograph then its complement G is also a cograph 

DEFINITION: A cotree is a tree representing a cograph. The leaves of the cotree represents 
the vertices of the cograph. The internal nodes are numbered either 0 or 1. The root is 
numbered 1 only. If a node is numbered 1 then all of its children are numbered 0 and 
vice-versa. There is an edge between two nodes if and only if their lowest common ancestor 
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is numbered 1. 

Cotree construction from the parse tree or definition tree is given in [28]. Every permu- 
tation graph can be represented by a permutation [25]. As the set of cographs is a proper 
subset of set of permutation graphs, a cograph can also be represented by a permutation. 

6.2 Permutation Representation 

A parallel algorithm for recognizing permutation graphs and constructing their permutation 
representation takes O(log 2 n) time with n 4 processors [29], where n is the number of nodes 
in the permutation graph. However, if cotree of a cograph is given then its permutation can 
be constructed in O(logn) time with n processors [28]. In this chapter, an optimal O(logn) 
time algorithm with processors is described. 

Let G = (V, E) be a cograph with cotree T. Let L v be the set of leaves in the subtree 
rooted at v. The algorithm of [28] plots the vertices V of G on an n x n grid. For each 
node v of T, the vertices of G v ( the subgraph induced by vertices of Z„), are plotted in 
| L v | X | L v | square region on the grid. The coordinates of the lower left corner (x„, y v ) 
and the upper right corner ( x ' v , y' ) can be calculated in an inductive manner. 

If coordinates of v are known and if Uj, vj, . . . , w* &rc children of v then the corresponding 
coordinates of its children are computed as follows [28]: 


Case 1 — v is a 0 node See Figure 6.1(a): 

i— 1 

*Vi = + £ | 

1 

(6.1) 

j= 1 

i-1 

y Vi = Vv + J2\ 

1 

(6.2) 

3=1 

t 

< = + J2 1 

1 1 -i 

(6.3) 

j=l 

i 

y'vi = v* + 1 

K, 1 -i 

(6.4) 




Case 2 — v is a 1 node See Figure 6.1(b): 

»-i 

Xvj = x v + ^2 | Zvj | 
3 = 1 


(6.5) 
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v« = sv + E 1 1 c I (#•«) 

j=t-i 

< = *» + E|£. i l- 1 (6-7) 

i - 1 

< = ». + E|£, i |-l (6-8) 

j=* 

For a leaf node to, the lower-left-corner and the upper-right-corner coincide; i.e. x w = x' w 
and y w = y' w . The permutation representation for the cograph, ir can be obtained as follows: 
for each leaf w let ir Xw = y w . 

The recursive equations 6.1, 6.2,. ..,6.8 can be solved in parallel as follows: 

1) For node v with children ...,«fc compute the prefix and the postfix 

summation of the corresponding | L Vi |’s. 

2) Carry out prefix computation on the cotree (see Chapter 4 for details). 

As both these steps can be implemented in O(logn) time with (j~) processors on an 
EREW PRAM, the following theorem holds: 

Theorem 6.1 Permutation representation of a cograph, given its cotree representation, can 
be obtained in O(logn) time with (j^j) processors on an EREW PRAM. I 

6.3 Coloring 

Algorithms for maximum weighted clique and minimum coloring of a permutation graph of 
[2]. take O(log 2 ) with n 3 processors. However, the maximum clique of a cograph (given its 
cotree) can be obtained in O(log n) time using processors on an Exclusive Read and 
Exclusive Write PRAM [1]. The algorithm given in [1] can be extended to find a maximum 
weighted clique; as cograph is a perfect graph, the chromatic number is equal to its clique 
number. We obtain an algorithm for minimum coloring with the same resourse bounds. 

The algorithm of [1] is first described for completeness; later the algorithm is extended 
to obtain a minimum coloring. Let G v be the subgraph induced by the vertices of L v . And 
let C(v) be the chromatic (or clique) number of G v . The chromatic number of each vertex 
in a cotree can be calculated as follows: 
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Case 1 — v is a leaf node: chromatic number of single vertex graph is one, hence C( v) 
is equal to one in this case. 

Case 2 — v is numbered 0: Let v has k children, namely t>i, v 2 , . . ., v*. As v is numbered 

0, G v is disjoint union of G Vl , G ^ , . . . , G Vk . Hence C(v) will be max(C(t>i ), C(v 2 ), ...,C(vk 

Case 3 — v is numbered 1: Let v has k children, namely vi, v 2 , . . ., ujt- As v is numbered 

1, in G v there will be all possible edges between subgraphs G Vl , G V2 , . . . , G Vk . Hence 
chromatic number of C(v) will be summation of chromatic number of its children. 

Theorem 6.2 ([1]) Chromatic number of a cograph can be found in 0(logn) time using 
processors on an EREW PRAM model. ® 

If C(v) is m then we need m colors to color the subgraph G v . Let us say we color G v 
with colors [i . . . C +m — 1]; we say that the range of node v is [i . . .* + m — 1]. Let the 
chromatic number of root, r be c, i.e. the given cograph is colorable with c colors. Let the 
root’s children be vi,v 2t . .., v* with ci,c 2 , . . . ,cjt as their chromatic numbers respectively. 

As root is numbered 1, we know that c = c\ + c 2 d — • + £*. The ranges of each child of root 
is computed as below: 

range of vi will be [1 . . . ci] 

range of v 2 be [1 + ci . . .c\ + c 2 ] 

range of u,- will be [1 + c i ' ' '53j=i c iL for 3 < i < k 

The above formulae are true for any node numbered 1 except that instead of 1 in the 
formulae it should be its lower value. Let the range of a node v be p . . .q. the ranges of 
children of v will be calculated as follows. 

range of vj will be [p . . . Ci] 

range of r 2 will be [p + Ci . . . C\ + c 2 ] 

range of w,- will be [p + 1 c i ' ’ '£j=i c il> f° r 3 < i < k 

If the node is numbered 0 then its range can be carried to all of its children. 
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For a loaf no<lo the tipper and lower values of the range coincide. In fact the upper 
value of any node can be p + C(v) — 1, where p is its lower value and C(v) is its chromatic 
number. Hence propagation of lower value is enough. 

Above procedure can be implemented efficiently as follows: 

1) Find C(v), the chromatic number of each v in the cotree. 

2) Carry out the following sub steps: 

a) Let for each node v numbered 0, /3(v) = Yii C(b{) — C(v) where 
bi is a right brother 1 of v. 

b) Let for each node v numbered 1, /3(v) := 0 

3)Apply prefix computation on trees (see Chapter 4) for /3(v) with *+’ as the 
associative binary operator. 

Hence the following theorem follows, 

Theorem 0.3 Minimum coloring of a cograph can be done in 0(logn) time using (j^j) 
processors on an EREW PRAM model if its cotree is given. I 


J If Vl ,t> 2 ,... , vh are 


children of a node v in that order, we say, vj, v a , . . . , v< we right brothers of node 



Chapter 7 


Conclusions and Scope for Further 
Work 


In this thesis efficient parallel algorithms have been obtained for some of the interesting prob- 
lems on Special graphs like Interval graph, Permutation graph, Cograph etc. It is possible 
to design better algorithms on these graphs because they possess some special properties 
which can be used to obain better algorithms. An interesting problem not considered in 
this thesis is that of recognition of these graphs. 

Test for connectivity of an interval graph, given its interval model, takes at least 
0(n log n) time where n is the number of intervals [39]. However, if the intervals were given 
in sorted order then it is possible to design even doubly-logarithmic time and even faster 
algorithms using k-th-prefix-maxima technique. It may be possible to use this technique to 
obtain better algorithms for other problems also. 

However, this technique does not work for k edge connected components. To the best 
of our knowledge there is no known algorithm for finding k edge connected components of 
an interval graph. Finding such an algorithm is an interesting open problem. 

The best algorithm (sequential) for finding maximum weighted independent set of a per- 
mutation graph takes O(nloglogn) time [52]. The fastest parallel algorithm takes 0(log 2 ) 
time [52] and is far from being optimal. The problems of finding optimal O(logn) time par- 
allel algorithms for maximum weighted independent set and minimum weighted dominating 
set are open. 
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Appendix A 


Strong Perfect Graph Conjecture 


A.l p-Critical Graphs 

An undirected graph G{V, E) is called p-critical if G(V,E ) is not perfect but every proper 
induced subgraph of G is perfect. In other words removal of any vertex makes the graph 
perfect. More formally, for any node v the following are true: 

u(G - v) = ~ ») 

a(G-v) = k(G-v) 

where G — v represents the subgraph induced by all vertices of V but v 
Theorem A.l (see [25]) // G(V,E) is a p-critical graph, then 

a(G)co(G) + 1 =| V | and for all vertices v of G, a(G ) = k(G — v) and co(G) = x(G — v). 

The odd cycle Cik+i x , where k > 2, is not a perfect graph since its clique number is 2 
and its chromatic number is 3. But removal of any node from C 2 fc+i results in a perfect 
graph or in other words every proper subgraph of Cik+i is perfect. Its complement C^fc+i 
is also p-critical. In fact, C 2 fc+i and Cik+i (where k > 2) arc the only p-critical graphs 
known to date. 

Berge conjectured that there are no other p-critical graphs. This is known as Stivng 
Perfect Graph Conjecture (SI’GC). 

‘Cjk+i is a connected graph of 2fc + 1 vertices in which every vertex degree is 2. 


(A.1) 

(A.2) 
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A. 2 Equivalent forms of SPGC 

SPGC-l: An undirected graph G is perfect iff it contains no induced subgraph isomorphic 
to Cjjfc+i or C^k+i where k > 2. 

SPGC-2: The only p-critical graphs that exist are Cjfc+i or C 2Jt+i where k>2. 
SPGC-3: there is no p-critical graph with a(<j) > 2 and oj(G) > 2. 

If we restrict the universe of graphs being considered, then in some cases the SPGC can 
be shown to hold. It is proved that the SPGC holds for the following problems: 

• Planar graphs [49]. 

• K i ,3- free graphs [4 0] . 

• Graphs whose chromatic number is less than or equal to 3 [50]. 

• Toroidal graphs; graphs having maximum vertex degree lee than or equal to 6 [26]. 



Bibliography 

[1] Abrahamson, Dadoun, Kirkpatrick and Przytycka. A simple parallel tree contraction 
algorithm. Journal of Algorithms, 10(l):287-302, 1989. 

[2] Adhar G. S. and Peng S. T. Parallel Algorithms for Cographs and Parity Graphs with 
Applications. Journal of Algorithms, ll(2):252-284, 1990. 

[3] Aho A., Hopcroft J. and Ullman J. The Design and Analysis of Computer Algorithms. 
Addison-Wesley, 1974. 

[4] Albers, S., Hagerup,T. Improved Parallel Integer Sorting without Concurrent Writing. 
In Proc 3rd Ann. ACM-SIAM Symp on Discrete Alg., pages 463-472, 1992. 

[5] Alt II., Hagerup T., Melhorn K., and Preparata F. Deterministic simulation of ide- 
alised parallel computers on more realistic ones. In proc. Mathematical Foundations of 
Computer Science, Lee. Notes in Comp. Science, Vol. 233. Springer- Verlag, 1986. 

[6] Anderson R. J. and Miller G. L. Deterministic Parallel List Ranking. Algorithmica, 
6(6):859-868, 1991. 

[7] Beame, P., Hastad, J. Optimal Bounds for Decision Problems on the CRCW PRAM. 
Journal of ACM, 36(3):643-670, 1989. 

[8] Berkman 0. Paradigms for Very Fast Parallel Algorithms, Pli. D. thesis, Tel Aviv 
University. UMIACS-TR-91-117, University of Maryland, August 1991. 

[9] Berkman O., Ja’Ja’ J., Krishnamurthy S., Thurimella R. and Vishkin U. Some Triply- 
Logarithmic Parallel Algorithms. In Proceedings IEEE FOCS, pages 871-881, 1990. 


43 



44 


{10] Berkman 0., Schiebcr B. and Vishkin U. Optimal Doubly Logarithmic Parallel Algo- 
rithms Based on Finding ail Nearest Smaller Values. Journal of Algorithms, 14(3):344- 
370, 1993. 

[11] Bertossi A. A. and Bonuccelli M. A. Some Parallel Algorithms on Interval Graphs. 
Discrete Applied Mathematics, 16:101-111, 1987. 

[12] Bhatt,P.C.P., Diks,K., Hagerup,T., Prasad, V.C., Radzik,T., Saxena,S. Improved De- 
terministic Parallel Integer Sorting. Information and Computation, 94:29-47, 1991. 

[13] Cherian J ., Kao M. Y. and Thurimella R. ‘Scan first Search and Sparse certificates: An 
improved parallel algorithm for k-vertex connectivity. SIAM Journal of Computing, 
22(1):157-174, 1993. 

[14] Chiang Y. and Fu K. S. Parallel processing for distance computation in syntactic 
pattern recognition. In Proceedings of IEEE Computer Society Workshop on Computer 
Architecture for Pattern Analysis and Image database management, 1981. 

[15] Cole R. and Vishkin U. Deterministic Coin Tossing with Applications to Optimal 
Parallel List Ranking. Information and Control, 70(l):32-53, 1986. 

[16] Cole R. and Vishkin U. Approximate parallel scheduling part I: the basic technique 
with applications to optimal parallel list ranking in logarithmic time. SIAM Journal 
of Computing, 17(1):128-142, 1988. 

[17] Cook S., Dwork C. and Reischnk R. Upper and lower time bounds for Parallel Random 
Access Machines without Simultaneous writes. SIAM Journal of Computing, 15(1):87- 
97, 1986. 

[18] I). G. Corneil, Y. Perl and L. K. Stewart. A Linear Recognition Algorithm for Cographs. 
SIAM Journal of Computing, 14(4):926-934, 1985. 

[19] Das S. K., Calvin and Chen C. Y. Efficient Parallel Algorithms on Interval Graphs. In 
PARLE 92, Lee. Notes in Comp. Science, Vol. 605, pages 131-141. Springer- Verlag, 
1992. 



45 


[20] D<x> N. Graph Theory with Applications to Engineering and Computer Science. Scries 
in Automatic Computation. Prentice-llall, 1974. 

[21] Even S., Pnueli A. and Lempel A. Transitive Orientation of graphs and identification 
of Permutation Graphs. Canadian Journal of Mathematics, 23:160-175, 1971. 

[22] Even S., Pnueli A. and Lempel A. Permutation graphs and Transitive graphs. Journal 
of ACM, 19(3):400— 410, 1972. 

[23] Faith I'., Prabhakar 11. and Wigderson A. Relations between concurrent-write models 
of computation. SIAM Journal of Computing, 17(3):606-627, 1988. 

[24] Garey M. R. and Johnson D. S. Computers and Intractability: A Guide to the Theory 
of NP-completeness. Freeman, San Francisco, 1979. 

[25] Golumbic M. C. Algorithmic Graph Theory and Perfect Graphs. Computer Science 
and Applied Mathematics. Academic Press, New York, 1980. 

[26] Grindstead C. M. Toroidal graphs and the strong perfect graph conjecture. PhD thesis, 
University of California, Los Angeles, 1978. 

[27] Harary F. Graph Theory . Addison- Wesley, 1968. 

[28] He X. Parallel Algorithm for Cograph Recognition with Applications. Journal of 
Algorithms , 15(2):284-313, 1993. 

[29] He! in bold I). and Mayr E. Perfect graphs and parallel algorithms. In Proceedings, 
International Conference on Parallel Processing, pages 853-860, 1986. 

[30] Jung II. A. On a Class of Posets and the Corresponding Comparability Graphs. Journal 
of Combinatorial Theory, Series B , 24(2):125-133, 1978. 

[31] Kim H. Finding maximum independent set in a permutation graph. Information 
. Processing Letters, 36(l):19-23, 1990. 

[32] Knuth D. E. The Art of Computer Programming, volume 1. Addison- Wesley, 1973. 



46 


[33] Kruskal C. P., Rudolph L. and Snir M. The power of Parallel Prefix. IEEE Transactions 
on Computers , C-34(10):965-968, 1985. 

[34] Liang Y., Dhall S. K. and Lakshmivarahan S. Efficient Parallel Algorithms for Finding 
Biconnected Components of Some Intersection Graphs. In ACM Computer Science 
Conference , pages 48-52, 1991. 

[35] Lovasz. A characterization of perfect graphs. Journal on Combinatorial Theory, Series 
B, 13:95-98, 1972. 

[36] Lovasz. Normal hyper graphs and the perfect graph conjecture. Discrete Mathematics , 
2:253-267, 1972. 

[37] Lu Mi. Parallel Computation of Longest-Common-Subsequence. In Advances in 
Computing and Information, Lee. Notes in Comp. Science, Vol. 468, pages 385-394. 
Springer- Verlag, 1990. 

[38] Modelev8ky J. L. Computer applications in applied genetic engineering. Advances in 
Microbiology, 30:169-195, 1984. 

[39] Olariu S., Schwing J. L. and Zhang J. Optimal Parallel Algorithms for Problems 
Modeled by a Family of Intervals. IEEE transactions on Parallel and Distributed 
Systems, 3(3):3f>4-374, 1992. 

[40] Parthasarathy K. It. and Ravindra G. The strong perfect graph conjecture is true for 
ki t 3 -free graphs. Journal of Combinatorial Tlleory, scries B, 21:212-223, 1976. 

[41] Prabhakar R. The parallel Simplicity of Compaction and Chaining. Journal of Algo- 
rithms, 14(3):37 1-380, 1993. 

[42] Ramkumar G. D. S. and Rangan C. P. Parallel Algorithms on Interval Graphs. In 
Proc. hit. Conf. Parallel process., pages 72-74, 1990. 

[43] Saxena S. Lecture Notes for course CS742, 1992-1993, Second Semester. Dept, of 
computer science and engineering, IIT Kanpur. 



47 


[44} Spinrad J. On Comparability and Permutation Graph. SIAM Journal of Computing, 
14(3):658-670, 1985. 

[45} Sprague A. and K ulkarni K. 11. Optimal parallel algorithms for finding cut vortices 
and bridges of interval graphs. Information Processing Letters, 42(4):229-234, 1992. 

[46} Springsteel and Stojinouovic. General parallel prefix compulation with geometric, 
algebraic, and other applications. International Journal of Parallel Programming, 
18(6):485-5Q3, 1989. 

[47} Sridhar M. A. and Goyal S. Efficient Parallel Computation of Hamiltonian Paths and 
Circuits in Interval Graphs. In Proc. Int. Conf Parallel process., pages 83-90, 1991. 

{48} Sumner D. P. Dacey graphs. J. Austral. Math. Soc ., 18(4):492-502, 1974. 

[49} Tucker A. C. The strong perfect graph conjecture for planar graphs. Canadian Journal 
of Mathematics, 25:103-114, 1973. 

[50} Tucker A. C. Critical perfect graphs and perfect 3-chromatic graphs. Journal of 
Combinatorial Theory, series B, 23:143-149, 1977. 

[51} Wagner H, Triangulating a monotone polygon in parallel. In Computational Geometry 
and Its Applications, Lee. Notes in Comp. Science, Vol. 383, pages 136-147. Springer- 
Verlag, 1988. 

[52} Yu M. S., Tseng L. Y. and Chang. Sequential and parallel algorithms for the maximum- 
weight independent set problem on permutation graphs. Information Processing Let- 
ters , 46(1):7-11, 1993. 



